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ABSTEIACT 


Some  of  the  results  contained  in  this  paper  have  been  obtained 
as  a  by-product  of  other  research,  but  have  been  collected  in  this  re¬ 
port  to  establish  a  frame  of  reference  part  of  which  is  applicable  to 
flame  pattern  configuration.  The  first  part  of  the  paper  is  directed 
towards  investigating  the  distribution  of  casualty  producing  material 
from  the  dynamics  of  dispersion.  Some  of  this  information  may  prove 
useful  in  establishing  theoretical  upper  bounds  on  the  region  of 
Lethality. 

The  last  part  of  the  paper  is  concerned  with  an  analytical 
model  for  flame  distributions  resulting  from  a  bomb  dispersed  gel, 
and  of  a  flamethrower  as  based  on  empirical  data.  This  data  will 
serve  as  a  lower  limit  of  desired  effects  since  it  is  currently 
obtainable.  A  number  of  curves  are  included  showing  the  area  of 
the  target  covered  as  a  function  of  the  orientation  of  the  strike 


relative  to  the  target  center. 


CRITERIA  FOR  CASUALTY  PRODUCTION 

1.  Statement  of  Problem: 

There  are  several  facets  of  casualty  production  that  must  oe  in¬ 
vestigated  if  a  sensible  evaluation  jf  criteria  is  to  be  reached.  The 
following  discussion  will  be  directed  towards  the  problem  of  casualty 
production  as  a  consequence  of  non-nuclear  weapons. 

The  first  facet  of  the  problem  is  to  establish  a  criterion  of 
casualty  production  through  the  determination  of  parameters  necessary 
for  a  particle  to  break  the  surface  of  the  akin  either  by  impact  and 
penetration,  or  by  burning  This  problem  will  involve  questions  of 
size,  shape,  and  possibly  velocity  of  the  particle  as  well  as  questions 
of  protective  shielding,  be  it  clothing  or  some  kind  of  armour. 

Closely  related  to  the  foregoing  is  the  associated  problem  of 
establishing  criteria  for  casualty  production  as  dependent  upon  the  type 
of  activity  involved  and  the  nature  or  location  of  the  wound.  Detailed 
investigation  of  this  problem  will  confront  the  investigators  with  psy¬ 
chological  characteristics  of  the  personnel  involved  as  related  to  the 
tactical  situation  as  well  as  physiological  criteria.  Further  consider¬ 
ation  of  this  problem  will  not  be  considered  here. 

A  third  problem  has  to  do  with  the  distribution  of  the  material 
effecting  the  casualty  production  about  the  point  of  impact,  that  is, 
the  establishment  of  a  Region  of  Lethality  and  the  associated  density 
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of  casualty  production  within  this  region.  This  then  relates  to  the 
second  problem  through  a  consideration  of  the  attitude  of  the  personnel. 

The  introductory  comments  of  this  discussion  will  be  rather 
general,  but  is  intended  as  a  springboard  for  outlining  more  specific 
requirements  without  producing  an  extensive  lacuna  or  gap  in  concepts. 

2 .  Parameters  of  Casualty  Production: 

As  a  working  definition,  we  shall  consider  an  individual  to  be  a 
casualty  if  he  is  unable  to  effectively  use  the  equipment  at  hand  due 
to  physical  injury  that  he  has  sustained.  The  most  tangible  form  of  in¬ 
jury  that  can  be  measured  physiologically  would  be  a  lesion  of  some  de¬ 
scription.  Lesions,  in  .urn,  can  be  produced  by  particles  penetrating 
the  skin  through  impact  or  by  burning. 

In  constructing  any  model  of  casualty  production  a  criterion  of 
mathematical  simplicity  is  desired.  However,  simplicity  to  permit  ease 
of  calculation  must  be  constrained  through  the  use  of  logical  connections 
of  the  employed  parameters. 

Any  criteria  of  casualty  production  should  be  stated  in  terms  of 
a  probability  distribution  that  will  increase  monoconically  in  terms  of 
of  the  critical  parameters. 

Ballistic  wounding  as  defined  by  the  piercing  or  breaking  of  the 
skin  through  impact  will  be  dependent  on  dynamical  parameters 

(1)  momentum,  mv 

2 

(2)  energy,  ^  mv 

(3)  density  of  particles  in  the  neighborhood  of  the  personnel 
By  considering  the  parameters  of  mass,  m,  and  velocity,  v,  in  their 


nascent  states,  illogical  connections  between  parameters  may  arise  which 
in  no  sense  can  be  theoretically  justified. 

The  following  conditions  can  be  expected  to  hold  in  any  probability 
function  for  casualty  production  by  "ballistic  wounding": 

(1)  If  the  momentum  and  energy  are  too  small,  there  will  be  no 
casualty  production  implying  a  threshold  value  below  which 
the  probability  of  wounding  is  zero.  In  other  words,  there 
is  a  lower  bound  ot  momentum  and  energy  below  .jhich  a  casu¬ 
alty  cannot  be  produced. 

(2)  As  momentum  and  energy  increase,  the  probability  of  wound¬ 
ing  will  increase  to  a  maximum  of  one,  provided  the  density 

of  wounding  particles  remains  sufficiently  high. 

,  2 

(3)  Since  momentum  M  «  mv  and  energy  E  ■  ^  mv  ,  any  increase 
in  one  of  the  parameters  implies  an  increase  in  the  other. 
Hence,  the  probability  should  be  monotone  increasing  with 

m  and  v . 

(4)  As  the  distance  from  the  center  of  impact  increases,  the 
density  of  wounding  particles  decreases  and  hence  the  prob¬ 
ability  of  wounding  decreases. 

In  considering  lesions  caused  by  "flame  smterials"  a  slightly 
different  framework  must  be  considered  in  that  the  dynamical  feature 
of  momentum  and  energy  of  oKition  are  not  of  maxisml  importance,  but 
density  of  flammable  material  is.  For  flame  materials,  it  is  required 
that 

(1)  size  of  "flame"  supporting  particle  is  large  enough  to 

burn  for  sufficient  tia*  at  lesion  producing  temperatures. 
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(2)  density  of  "flame"  supj  ’‘ting  material  be  sufficient  so  that 
fire  prevention  and  control  techniques  are  reduced  in  effec¬ 
tiveness,  preferably  to  a  point  of  being  completely  ineffec¬ 
tive  . 


3.  Probability  Function  Forms  for  Casualty  Production: 

From  the  foregoing  discussion  of  parameters  it  becomes  evident 
that  the  probability  function  of  "ballistic  wounding"  and  "flame  lesions" 
as  criteria  for  casualty  production  are  fundamentally  different. 

Since  "ballistic  wounding"  depends  upon  the  dynamical  parameters 
of  momentum  and  energy,  it  is  conjectured  that  a  probability  density  de¬ 
pendent  only  on  them  should  assume  the  form  of 

2  2 

a(mv  -  mv  )  +  2b(^  mv  -  ^  mv  ) 

P  *  1  -  e  °  ° 

or  more  simply 

mja(v  -  V  )  +  b(v  -  v  ^)  ) 

P  «  ’  °  °  - 

r  *  A  -  e 

Where  a  is  the  coefficient  of  momentum  change  and  2b  is  the  coeffi¬ 
cient  of  energy  change  corresponding  to  a  particle  of  a  given  mass  m. 

The  velocity  is  a  threshold  velocity,  below  which  no  wounding  takes 

place  for  that  size  of  particle.  If  sufficient  data  is  available,  the 
coefficients  a  and  b  can  be  estimated  from  a  least  squares  criterion 
applied  to  the  equivalent  linear  equation 

2 

ln(l  -  P)  •  ai^v  +  bm^v 

2  2  2 

Av«v-v,21v  •v-v 


where 
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If  the  dependence  on  distance  from  the  center  of  impact  is  to  be 

incorporated  in  the  probability  statement,  it  will  have  to  be  multiplied 

2 

by  some  factor  near  Gaussian  in  form  such  as  exp(-r  ) ,  The  true  nature 
of  this  factor  will  depend  upon  the  dynamical  paths  of  the  particles. 

The  probability  function  for  causing  "flame  lesions"  is  funda¬ 
mentally  a  time  and  density  dependent  function  and  hence  P  =  P(d,  t) . 

2 

For  example  it  has  been  determined  chat  20  gms/m  of  white  phosphorus 
are  needed  to  produce  casualties. 

4 .  Region  of  Lethality 

The  parameters  of  concern  for  the  production  of  casualties  are 
different  for  ballistic  wounding  and  wounding  caused  by  flame  lesions. 

The  determination  of  the  region  of  lethality  resulting  from  the  burst 
of  a  conventional  H.E.  warhead  can  be  approximated  from  certain  theo¬ 
retical  concepts.  In  order  to  obtain  a  solution  in  closed  form,  however, 
certain  s  mplifying  assumptions  need  to  be  made.  Some  of  these  same 
assumptions  are  plausable  for  determining  the  distribution  of  fuel  from 
a  flame  projectile.  Possible  modeling  structures  for  flame  warheads 
could  have  a  controlled  H.E.  burst  to  distribute  the  flaae  fuel  in  a 
prescribed  pattern  from  an  aerial  burst  or  a  projectile.  A  hybrid  pro¬ 
jectile  would  also  have  the  goal  of  the  K-  E.  fracturing  or  otherwise 
disrupting  the  target  surface  enabling  the  fuel  to  more  effectively 
ignite  the  target  material. 

It  is  assumed  that  the  particles  of  the  burst  are  finite  in 


number  (N) ,  have  constant  mass  (m) ,  and  that  their  density  or  dis¬ 
persion  will  reach  a  lower  limit  of  effectiveness.  It  is  furthermore 
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assumed  that  the  particle  distribution  will  be  governed  by  gravitational 
attraction  and  air  resistance.  The  altitude  of  the  burst  will  be  con¬ 
sidered,  but  assumed  to  be  near  enough  to  ground  level  so  that  the  ac¬ 
celeration  of  gravity  will  be  effectively  constant. 

The  differential  equation  for  the  motion  of  a  typical  particle 
in  vector  form  is 

R  =  F(R)  +  g£^, 

where  superposed  dots  denote  time  derivatives  and  capital  letters  denote 
vectors.  In  particular  R  is  the  position  vector,  is  a  unit  vector 

directed  downward,  g  the  acceleration  of  gravity  assumed  to  be  constant 
and  F  a  vector  function  of  the  velocity.  The  true  nature  of  the  vector 
F  is  not  known  explicitly,  however  its  behavior  lies  somewhere  between  a 
linear  response  to  velocity  and  a  quadratic  response  to  velocity. 

5 .  Air  Resistance  Proportional  to  the  Velocity; 

If  the  vector  function  F  is  proportional  to  the  velocity  the 
equation  of  motion  becomes 

R  =  -  kR  +  gEj^ 

This  has  the  first  integral  of 

F  -  -  [(gEj,  -  -  gsj  /k 

where  is  the  initial  velocity  of  the  particular  particle  or  glo¬ 

bule  when  the  warhead  explodes.  This  equation  yields  the  second  inte¬ 
gral  of 

R  -  [(gE^  -  kVpCe"*'*'  -  1)  +  k(gt  -  hk)Ej  IV? 
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where  h  is  the  height  above  the  ground  of  the  burst.  That  is  at  t  =  0 
R  .  Ej  -  -h. 

5a.  Region  of  Lethality  for  Ballistic  Wounding; 

To  inflict  damage  to  the  target  by  ballistic  wounding,  it  is  as¬ 
sumed  that  the  magnitude  of  velocity  of  a  particle  must  exceed  a  given 
velocity.  In  other  words,  if  IrI  <  v^ ,  there  will  be  no  effective 
damage  to  the  target.  The  lethal  region  is  bounded  by  the  envelope  of 
the  set  of  points  for  which  Ir!  =  v^ ,  provided  of  course  that  the  parti- 
cle  density  is  sufficiently  high. 

Each  particle  in  the  burst  will  define  a  separate  trajectory  which 
will  be  determined  by  its  own  initial  condition,  i.e.,  the  vector  velocity 
of  the  warhead  and  the  direction  and  velocity  that  the  particle  takes 
relative  to  the  warhead.  The  geometry  and  distribution  of  material  on 
the  warhead  will  affect  the  region  through  density  considerations. 

On  the  boundary  of  the  Lethal  region  R  »  Vj^,  from  which  it 
follows  at  the  boundary  that 

''l  +  =  ''o  * 

where  t^  is  the  time  for  the  particular  particle  to  reach  the  boundary 
and  Rg  is  the  position  vector  denoting  the  boundary.  It  also  follows 
that 


gE  -  kV^  -ktj  gEj 
e 


^-k 


"L  k 

Eliminating  V.  from  equations  (1)  and  (2)  gives  a  vector  equation 

id 

relating  R^  and  t^ 


(2) 
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kR 


B 


gE,  -kt 

(Vq  -  -^)(1  -  e  +  (gtg  -  hk) 


(3) 


The  initial  velocity  V^,  consists  of  two  parts 
i)  the  velocity  contributed  by  the  projectile 


V  =  V  (E,  sin  (I)  +  Eo  cos  uj) 

p  p  1  2 


(4) 


li)  and  the  Initial  velocity  imparted  to  the  individual  particle  on  ex¬ 
ploding  the  warhead 


*  v^  (Ej^  cos  a  +  E2  cos  B  +  cos  y)  (5) 

where  Ej^  is  a  unit  vector  oriented  downward,  E2  a  unit  vector  in  the 
horizontal  direction  of  the  direction  of  motion  of  the  projectile  and 
E^  a  unit  vector  perpendicular  to  Ej^  and  E2>  (Ej^  x  E2  *  E^)  . 

Let  the  vector  equation  of  the  boundary  be 


Rj  .  aEj  +  rE^  +  dEj 


(6) 


where  a  is  the  altitude,  r  the  range  and  d  the  lateral  displace¬ 
ment  of  the  particle  from  the  point  of  explosion.  Equation  (3)  and 
(6)  now  gives 


ka 


(v  sin  oj  +  V  cos  a 


-  h 


(1  - 


-ktB 

e  ®)+gt3 


hk 


-kt. 


kr 


(v  cos  U) 
p 


v^  cos  3)  (1  -  e 


(7) 


-ktB 

kd  ■  v^  cos  Y  (1  -  e  ) 

The  proportionality  constant  k  can  be  eliminated  from  two  of  these 
equations  as  well  as  the  exponential  term  to  yield 
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rv  cos  Y  ®  d(v  cos  uj  +  v  cos  B) 
e  '  p  e 


(8) 


Another  relation  independent  of  the  exponential  term  is 

(ka  -  gtg  +  hk)  cos  y  ®  ^ 

Solving  this  equation  for  t  yields 

D 


'B 


,  ,  ,  , .  j  d(v  sin  u)  +  V  cos  a) 

k(a  +  h)  ^  d _  ^  p  e _ _ 

g  kv^  cos  Y  gv^  cos  Y 


(10) 


I  /  L  u\  dC*"  -  (v  sin  u)  +  V  cos  &)  /  ) 
_  k(a  +  h)  ^  r  ^  p  e _  g 

g 


V  cos  V 
e  ’ 


k(.H-h>  ,  Ik  -  (V,  «)/J 

”  g  V  cos  (1)  +  V  cos  8 

B  p  *  e  ^ 


g  gk  (Vp  cos  u)  +  cos  p) 


As  a  first  approximation  1  -  e  =  kt  if  kt  <  1, 

D  D 

from  this  it  follows  that 


kd  ■  cos  Y  (1  ■  s  )  V*" 


B 


hence 


L  g 


d(v  sin  (I) 


kd  ■  kv  cos  Vi  <1 

e  ’  L  g  kv  cosy 


gVg  cos  Y 


+  V  cos  Of)! 

- - - I 

IS  V  J 


(11) 


It  now  follows  that 


k  (a  +  h)  v^  cos  y  ■  sin  «)  +  v^  cos  ») 


(12) 


Using  equation  (8)  and  (12)  in  the  identity 


2  2  2 
cos  Qt  +  cos  B  +  cos  Y  ■  1 


(13) 
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we  get 


dv  cos  (x  “  k(a  +  h)v'  cos  v  "  dv  sin  tu 
«  e  ^  p  * 

dv^  cos  B  =  rv^  cos  y  -  dv^  cos  u) 

Hence 


2  2  2  r  2  2  2*V  2  7 

d  (1  -  cos  y)  ■  l^k  (a  +  h)  j  Y  “  2dv^VpCOsy 

•  ^k(a  +  h)  sin  u)  +  r  cos  +  d^v  ^ 


This  can  be  rewritten  as 


2  2 
V  cos  Y 


+  r 


2dv  V  cosv 
e  p  ^ 

2  21 
'  Ve  J  . 


[k(a  +  h) 
■  0 


sin  u)  +  r  cos 


or 


0  V  cos  V 

e  ” 

\v  V  cos  y  sin  u)  -v  v  cos  v  cos  tu 
e  p  T  e  p  ^ 


(a  +  h 

r,  di 


■W 


This  is  a  homogeneous  quadratic  form  and  thus  represents  a  degenerate  conic. 
However,  an  assignment  to  the  initial  attitude  of  the  projectile  Assignment 
of  tu)  and  a  specification  of  y  which  determine  points  on  a  cone  with  its 
axis  horizontally  sideways,  constraint  on  the  boundary  values  of  a,  r,  and 
d  are  then  given.  By  further  assigning  a  value  to  d  compatible  with  these 
constraints  a  level  line  mapping  is  obtainable. 

Taking  d  as  a  parameter,  the  equation  can  be  rewritten  as 
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dv  sin  (1) 

2 

(a  +  h)  -  ,  ^ — 

'  '  kv  cos  V 

L  e  ’  J 

I 

^22,2 
d  V  sin  V 
e  ' 


,2  2  2 
k  V  cos  V 
e  ^ 

which  is  an  ellipse  with  its  center  at 


dv  cos 
^ _ 


V  cos  V 
e  ’ 


2  2  2 
d  V  sin  Y 
e  ^ 


V  cos  V 
e  ^ 


2 


(a  +  h,  r) 


^dVp  sin  0#  dVp  cos 
kv^  cos  Y  ’  cos  Y  j 


and  semi' axes  of  lengths 

d  tan  y  and  d  tan  y/^ 

in  the  r  and  a  +  h  directions  respectively. 
From  equatipn  (12)  find 


d(v  sin  (»  +  V  cos  or) 
''e  V  -  . ®kuT"h)  “ - 


which  is  then  substituted  into  the  equation  for  the  degenerate  conic 

|2 


r  - 


k(a  +  h)  ^  cos  tt) 


V  sin  u)  +  V  cos  Of 
P  ^ 


+  d 


f 

•  • 

V  sin  « 

f 

\f  .  .  n2 

~  V  sin  (1)  -f  V  cos  01 

► 

|(v  sin  u)  +  V  cos  or) 

p  e 

L  p  e 

d 

2  2 
V  cos  or 
e  _ 


j(Vp  sin  fl#  +  cos  a)‘ 


(v  sin  (t)  +  V  cos  or) 
'  p  e 


2  2 
k^(a  +  h)^ 


2  2  2  2 
k  V  ^  sin^  «  (a  +  h)^ 
e 


(v  sin  tt)  V  cos  or)' 
'  p  e 
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An  interpretation  of  this  equation  is  that  any  particle  emitted 
from  the  projectile  along  the  right  circular  cone  a  (axis  vertical) 
will  have  its  boundary  for  the  region  of  lethality  on  a  circle  that  has 
been  translated  forward  from  the  point  of  explosion.  For  a  given  h  and 
(U,  the  maximum  translation  occurs  with  a  -  TT/2.  The  translation  is  then 

k(a  +  h)  cos  0}  0)  0 

and  the  radius  of  the  boundary  circle  is 


kv  (a  +  h) 
e 


V  sin  (I) 
P 


U)  0 


The  ratio  of  the  translation  to  the  radius  of  dispersion  is  given  by 
Vp  cos  tu/v^  sin  a.  In  order  that  the  region  of  lethality  be  essentially 
circular  it  is  necessary  that  this  ratio  be  less  than  1  for  some  a.  In 


other  words  it  is  sufficient  if  v  cos  tu  <  v  to  have  a  circular  region 

P  6 

of  lethality. 


5b.  Region  of  Lethality  for  Flame  Lesions; 

In  order  to  produce  a  casualty  by  flames  the  velocity  of  the 
globules  of  burning  material  is  not  of  primary  concern.  Instead  we  are 
concerned  with  the  geometrical  distribution  of  the  material. 

Again  assuming  that  the  air  resistance  is  proportional  to  the 

*  •  • 

velocity,  the  foregoing  equations  for  R,  R,  and  R  still  apply.  We  are 
now  concerned  with  the  when  and  where  the  globules  of  material  strike 
the  ground.  The  material  will  reach  the  ground  when  R  •  ■  0.  Con¬ 

sider  for  simplicity  the  material  leaving  the  warhead  with  a  vector 
velocity 

•  aEj  +  b(cos  OB^  *  sin  OB^) . 
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Due  to  horizontal  synmetry  we  can  for  simplicity  assume  that  9  ■  0.  The 
material  will  now  reach  the  ground  when 

(g  -  ka)(e"^''-l)/k+  (gt-hk)  -  0 

The  distance  from  the  center  of  explosion  will  be  given  by  R  *  E2  evalu¬ 
ated  at  the  time  value  obtained  from  the  foregoing  equation  and  is  given 
by 

R.E2  -  -b(e‘^*'-l)/k  -  r 

Note  that  the  limiting  value  for  r  is  b/k. 

The  time  can  be  obtained  from  these  equations  in  terms  of  a,  b,  h  and  r 
to  get 

t  -  jgr  +  k(bh  -  ar^  /bg. 

Substituting  this  value  of  t  into  the  equation  for  r,  gives  on  re¬ 
arranging  terms 


Since  r  <  b/k  let  r  •  ^  number  n  >  1.  Then  a  value 

of  t  can  be  obtained  from 


in  particular  t.  ■ 


i 


r)- 


Substituting  the  assumed  value  for  r  and  the  calculated  t  Into  the 
time  equation 


tbg 

gives  the  value  of 


h 


kbh  ■¥  r(g  -  ak) 

needed  under  the  prescribed  conditions,  i.e. 
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In  firing  from  a  static  platform  a  «  0  we  obtain  the  following 

results 


n 

Ln  ■  kt 

T 

n 

1 

n 

k^h 

Ln  — 
n 

-  1  ■ 

2 

0.69315 

Ln  2 

1 

"  2  “ 

0.19315 

6.18080  ft/sec^ 

3 

0.40547 

Ln| 

1 

•  ~  as 

3 

0.07213 

2.30816  ft/sec^ 

5 

0.22314 

1  5 
Ln  ^ 

1 

H 

0.02. x4 

0.74048  ft/sec^ 

If  suitable  values  of  k  are  known  approximately  satisfying  the 
assumptions,  the  height  of  the  burst  required  for  distributing  the  flames 
over  prescribed  circles  can  be  approximated. 

6.  Air  Resistance  Proportional  to  the  Square  of  the  Velocity; 

Under  the  assumption  of  air  resistance  being  proporcionv'l  to  the 
velocity  squared,  the  differential  equation  can  be  written  as 

R  +  kvR  •  gEj^ 
or 

••  •  * 

R  kl  R  .  R  -  gE^ 

Where  T  is  the  unit  target  multiplying  by  Ri  gives 

#  •  %  • 

R*  R  -f  kR-T  R*R  - 


d  R«  B  •  *  • 

^  +  kR*T  R*R  -  gR*E 


which  in  turn  can  be  rewritten  as 


r  -2k/R*Tdt} 

Tt  2‘>'*'=  J 


-2k  jR'Tdt 


g  R'E  e 


/R*Tdt  «  /vdt  dt  •J^ 


ds  “  s 


hence  the  differential  equation  of  motion  becomes 


d  1  2 

—  —  V  e 

dt  2 


This  has  a  first  Integral 


.  -2ks' 

R'Ee 


. .  -2ks 
R'E  e  dt 


12  r  f-  1 

“e  p‘*'8y  E'E  ®  dtj 

The  parameter  k,  according  to  John  E.  Younger  (Advanced  Dynamics 
-  Ronald  Press,  1958,  p.  97),  Is  of  the  order  of  magnitude 

k  lO'^/ft. 

Some  further  approximations  can  be  made  towards  getting  the  ve¬ 
locity  explicitly.  On  rewriting  and  Integrating  by  parts 


f*  f  (  -2k8 

y  R.E  e  dt  •  ^  e  E*dR  •  J  e  E.Tds 

_  _  -2ks  ,  /•  -2ks 

-  It*  ♦hA' 


E»T 

-2ks 

2k 

6 

E-T 

-2ks 

2k 

6 

kE*N 


~ e  E*  |V'N  +  hN^ 


4k‘  4k‘ '  ‘ 

where  n  is  the  unit  normal  vector  and  primes  denote  differentiation  with 


respect  to  arc  length. 

If  the  motion  is  planar  N* 


-  hT  (N' 


tB  -  kN) 
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Furthermore  if  the  curvature  k  Is  approximately  constant 


At  t  ■  0,  s  ■  0,  V  •  v^;  T  and  N  must  be  specified  but  these 
depend  upon  the  dynamical  paths  of  the  material.  Ignoring  this  we  find 


This  particular  equation  of  motion  readily  yields  an  expression  for  the 
velocity  and  hence  could  be  used  for  obtaining  some  insight  in  "ballistic 
tiounding"  criteria.  However,  it  does  not  yield  a  vector  solution  for  the 
displacements  in  closed  form  and  consequently  tne  dimensions  of  the  distrl*' 
bution  a*  desired  for  flame  investigations  are  unknown. 
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In  order  to  use  the  foregoing  approximations,  explicit  statements 
for  R  would  be  necessary  in  order  to  get  appropriate  values  for  th<> 
curvature  k . 

7.  Special  Empirical  Configurations: 

Field  studies  have  been  carried  out  studying  the  effect  of  flames 
from  flame  throwers  and  Napalm  bombs.  The  primary  objective  of  one  oi 
these  studies  was  to  determine  measures  for  defense  against  conventiorial 
flame  weapons.  A  by-product  was  the  establishment  of  flame  patterns  re¬ 
sulting  from  them.  In  both  situations  it  was  observed  that 

(1)  A  flame  distribution  divides  into  a  central  high 
heat  zone  and  an  outer  splatter  zone. 

(2)  Convective  currents  funnel  the  flame  towards  the  center. 
Relative  dimensions  of  the  high  heat  zone  and  splatter  zone  have  been 
determined.  In  order  to  more  coiq>letely  assess  the  lethal  area  produced 
by  these  two  conventional  modes  of  flame  delivery  it  becomes  necessary 

to  assess  the  relative  areas  of  the  two  zones  under  a  probability  distri¬ 
bution.  In  a  technical  report  for  Weapons  Systems  Analysis  Branch,  USA 
CDC  Artillery  Agency,  Fort  Sill,  OklahomSfa  nomograph  was  prepared  for 
determining  the  common  area  between  two  overlapping  circles  distributed 
with  a  Gaussian  Distribution  about  a  given  aim  point.  As  a  first  order 
model  it  is  conjectured  that  the  flame  distributions  of  flasm  thrower 
and  Napalm  bombs  can  be  approximated  by  two  eccentric  circles.  The 
smaller  inner  circle  representing  the  high  intensity  zone. 

Two  separate  sndels  have  been  considered,  based  upon  data  from 
the  report  "Individual  Defense  against  Napalmf*.  In  this  report  it  was 
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observed  that  a  110  gal.  Napalm  bomb  distributed  in  a  high  intensity  zone 
nearly  elliptical  with  dimensions  of  10  by  65  meters.  The  boundary  of 
the  splatter  zone  was  contiguous  at  a  point  of  inq>act  and  had  dimensions 
of  30  by  100  meters  as  in  Figure  1. 


Figure  1 

In  the  same  report;  it  was  observed  that  a  flame  thrower  at  69 
feet  had  a  much  more  eccentric  pattern  in  the  high  intensity  zone  of 
4  ft.  by  30  ft.  surrounded  by  a  splatter  zone  of  15  ft.  by  90  ft.  This 
time,  however,  the  zones  were  more  nearly  concentric  as  in  Figure  2» 


In  order  to  model  these  configurations  we  make  the  following 


assumptions : 


19 

(1)  The  target  under  fire  Is  circular. 

(2)  The  flame  pattern  can  be  reasonably  approximated  by  suit* 
ably  oriented  circles. 

At  this  point  it  becooies  necessary  to  specinTize  to  the  different  patterns. 
7a.  Bomb  Configuration: 

In  constructing  this  model  it  is  again  observed  that 

A.  High  intensity  and  splatter  zone  are  contiguous  on  the  leading 
edge. 

B.  The  zone  areas  have  the  ratio 

High  Intensity:  Splatter  Zone  »  1:5 

from  which  it  follows  that  if  the  high  intensity  zone  has  a 
radius  of  R,  then  the  splatter  zone  has  a  radius  of  VT  R. 

G.  Normalize  the  configuration  by  assigning  the  target  area  to 
a  radius  of  one  (1). 

D.  The  burst  studied  will  be  distributed  with  a  circular  Gaussian 
distribution  about  the  target  center  for  the  high  intensity 
zone,  but  with  a  constant  direction  of  delivery.  See  Figure 
3. 

The  area  damage  effected  by  the  high  intensity  zone  depends  upon 
the  distance  from  the  center  of  the  target.  This  data  is  partially  docu¬ 
mented  in  Table  I  and  exhibited  in  Figure  5.  The  area  of  the  splatter 
zone  depends  on  the  relative  direction  from  the  target  center  as  well  as 
the  distance.  Agalns  referring  to  Figure  3,  one  observes  that  the  dis¬ 
tance  of  the  center  of  the  circle  representing  the  splatter  zone  is  given  by 
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Figure  3 


The  distances  (1-b)  were  obtained  by  interpretation  from  the  data  of 
Table  II  which  is  a  listing  of  values  obtained  from 

1  -  b  “ 

where  R  was  assigned  the  values  of  0,1,  0.3,  and  0.5.  The  angle  cp 
was  assigned  values  at  18°  Increments  ranging  from  0°  to  180°  and  a 
the  distance  from  the  target  boundary  was  assigned  increments  of  0.1 
ranging  from  zero  the  boundary  to  1,  the  target  center. 

For  a  bomb  type  configuration  the  relative  dimensions  of  the 
parts  to  Figure  3  are  given  by  Table  II. 

The  distance  from  the  center  of  the  target  j  to  the  center  of  the 
splatter  zone  in  terms  of  a,  the  distance  from  the  target  boundary,  r. 
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the  distance  between  centers  of  the  high  intensity  zone  and  splatter  zone 
and  the  relative  direction  from  the  line  of  fire  to  impact  center  is  given 
for  selected  value  of  r  in  Table  III. 

Oi ce  Lhe  distance  from  the  center  of  the  target  is  obtained,  the 
relative  area  of  the  target  covered  by  the  composite  pattern  can  be 
readily  obtained  from  Figure  5.  Area  coverage  curves  for  this  pattern 
are  exhibited  in  Figures  6-10  inclusive. 

In  order  to  obtain  the  lethality  of  the  composite  pattern,  the 
total  area  must  be  separated  into  two  distinct  parts.  To  the  high  in¬ 
tensity  zone,  a  high  probability  would  be  assigned,  but  the  splatter 
zone  would  have  a  small  probability.  It  follows  that  the  probability  of 
casualty  production  by  such  a  pattern  would  then  be 

P  *  *  Ps*s 

where  p^^  »  probabiitiy  of  a  casualty  in  the  high  intensity  zone 

p  »  probability  of  a  casualty  in  the  splatter  zone 
s 

*  area  of  high  intensity  zone  within  the  target 
A^  »  area  of  splatter  zone  within  the  target. 

It  is  conjf''’tured  that  from  this  data  a  modeling  pattern  can  be  derived 
for  investigating  the  composite  effects  of  several.  palm  type  bombs. 

7b.  Flame  Thrower  Configurations ; 

The  configurations  resulting  from  the  use  of  flame  throwers  depend 
to  some  extent  upon  the  manner  of  employment.  Since  the  target  is  ob¬ 
servable,  except  for  the  obscuring  resulting  from  the  flame  and  ensuing 
smoke,  the  effect  can  be  monitored  by  the  operator.  In  this  section  we 
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will  not  go  Into  much  detail  since  it  is  not  felt  to  be  as  necessary. 

The  primary  observed  difference  is  that  the  high  intensity  zone 
is  no  longer  contiguous  to  the  splatter  zone  boundary,  but  does  have  its 
boundary  at  the  center  of  the  splatter  zone. 

A  second  variation  lies  in  the  observed  ratio  of  areas,  l.e.. 
High  Intensity  :  Splatter  Zone  *  1  :  11 
Utilizing  the  same  simplifying  assumptions,  the  configuration  is 
typified  by  Figure  4.  The  relative  dimensions  in  Figure  4  are  given  in 
Table  III. 


Figure  4 
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Curves  exhibiting  the  composite  areas  of  coverage  are  given 
in  Figures  11  to  15  inclusive. 

8.  Summary  and  Recommendations: 

The  Lethality  patterns  as  described  in  some  detail  give  some 
indication  of  the  probability  of  casualty  production  for  a  single  burst, 
be  it  a  bomb  o^.  a  flame  thrower.  These  configurations  will  be  used  as 
a  guide  in  developing  pattern  effects  and  attempting  to  prescribe  ways 
of  optimizing  these  effects. 

The  data  of  this  analysis  will  be  subsequently  related  to  the 
distribution  function  resulting  from  the  delivery  mode  and  expected 
values  for  casualty  production  will  be  estimated. 

Subsequently,  flame  patterns  will  be  designed  to  give  reasonable 
approximations  to  the  composite  pattern  of  a  single  burst.  This  will 
be  analyzed  on  the  OURI  SADI  Mark  IV  to  determine  effects  resulting  from 
a  number  of  bursts  deployed  against  a  target.  The  SADI  Mark  IV  will  in 
particular  be  used  to  determine  statistical  limits  of  effectiveness  that 
can  be  anticipated  to  hold  when  delivery  errors  of  different  types  are 
hypothesized. 


Fractional  Part  of  Area  Destroyed 
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Figure  5  -  Distance  Between  Centers  of  Circles 
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High  Intensity  Zone  Radius  0.3 


Distance  Between  Centers  Splatter  Zone  Radius 

0.371  0.671 

Figure  8 


6 


DUtanct  Betw««n  C«nt«rs 
0.494 


Figure  9 


Splatter  Zone  Radius 
0.894 
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High  Intensity  Zone  Radius  0.5 


Distance  Between  Centers  Splatter  Zone  Radius 

0.618  1.118 


Figure  10 
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High  Intensity  Zone  Radius  0.1 


18  36''  54''  72''  90''  lOS"  126’  144’  162"  180' 


Distance  Between  Centers 

0.2 


''platter  Zone  Radius 
0.663 


Figure  12 


30 


High  Intensity  Zone  Radius  0.3 


Splatter  Zone  Radius 
0.995 


Distance  Between  Centers 
0.3 


Figure  13 
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High  Intensity  Zone  Radius  0.4 


Distance  Between  Centers  Splatter  Zone  Radius 

0.4  1.327 


Figure  14 


High  Intensity  Zone  Radius  0.5 


Distance  Between  Centers  Splatter  Zone  Radius 

0.5  1.659 

Figure  15 


TABLE  I 


Area  of  Overlap  Between  Two  Circles 
A(R,B) 


b\L 

0.1 

0.2 

0.3 

0.4 

0.5 

0.5 

0.01 

0.04 

0.09 

0.16 

0.2500 

0.6 

0.01 

0.04 

0.09 

0.1600 

0.233 

0.7 

0.01 

0.04 

0.0900 

0.1466 

0.206 

0.8 

0.01 

0.0400 

0.0788 

0.  ’29 

0.175 

0.9 

0.01000 

0.0316 

0.0612 

0.0983 

0.143 

0.95 

0.00802 

1.0 

0.00490 

0.0192 

0.0421 

0.0731 

0.112 

1.05 

0.00183 

1.1 

0.00000 

0.0073 

0.0240 

0.0491 

0.082 

1.2 

0,00000 

0.0000 

0.0088 

0.0275 

0.054 

1.3 

0,00000 

0.0000 

0.0000 

0.0099 

0.030 

1.4 

0.00000 

0.0000 

0.0000 

0.0000 

0.011 

1.5 

0.00000 

0.0000 

0.0000 

0.0000 

0.0000 

txA(R,B)  =  cos 


■1  B^  +  (1  - 


2B 


+  R  cos 


-1  B  -  (1  -  R  ) 
2BR 


-  \  +  R^)  -  p  -  (R  -  l^J 
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TABLE  II 


Bomb  Configuration  (See  Figure  3) 


Radius  of 

Distance  Between 

Radius  of 

High  Intensity 

Center  of  High  Intensity 

Splatter 

Zone 

and  Splatter  Zone 

Zone 

0.1 

0.124 

0.224 

0.2 

0.247 

0.447 

0.3 

0.371 

0.671 

0.4 

0.494 

0.894 

0.5 

0.618 

1.118 

TABLE  III 

Flame  Thrower  Configuration  (See  Figure  4) 


Radius  of 

Distance  Between 

Radius  of 

High  Intensity 

Center  of  High  Intensity 

Splatter 

Zone 

and  Splatter  Zone 

Zone 

0.1 

0.1 

0.332 

0.2 

0.2 

0.663 

0.3 

0.3 

0.995 

0.4 

0.4 

1.327 

0.5 

0.5 

1.659 
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Table  IV 


Distance  to  Center  of  Splatter  Zone 


0° 

18° 

36° 

54° 

72° 

90° 

108° 

126° 

144° 

162° 

180° 

a 

-  0.0 

r 

»  0.1 

1.10 

1.10 

1,08 

1.06 

1.03 

1.00 

0.97 

0.94 

0.92 

0.91 

0.90 

0.3 

1.30 

1.29 

1.26 

1.20 

1.13 

1.04 

0.95 

0.86 

0.78 

0.72 

0.70 

0.5 

1.50 

1.48 

1.43 

1.36 

1.25 

1.12 

0.97 

0.81 

0.66 

0.55 

0.50 

a 

=  0.1 

r 

=  0.1 

1.00 

1.00 

0.98 

0.96 

0.94 

0.91 

0.87 

0.85 

0.82 

0.81 

0.80 

0.3 

1.20 

1.19 

1.16 

1.10 

1.03 

0.95 

0.86 

0.76 

0.68 

0.62 

0.60 

0.5 

1.40 

1.38 

1.34 

1.26 

1.16 

1.03 

0.88 

0.73 

0.58 

0.45 

0.40 

a 

»  0.2 

r 

«  0.1 

0.90 

0.90 

0.88 

0.86 

0.84 

0.81 

0.77 

0.75 

0.72 

0.71 

0.70 

0.3 

1.10 

1.09 

1.06 

1.01 

0.94 

0.85 

0.76 

0.67 

0.58 

0.52 

0.50 

0.5 

1.30 

1.28 

1.24 

1.17 

1.07 

0.94 

0.80 

0.65 

0.49 

0.36 

0.30 

a 

«  0.3 

r 

-  0.1 

0.80 

0.80 

0.78 

0.76 

0.74 

0.71 

0.68 

0.65 

0.62 

0.61 

0.60 

0.3 

1.00 

0.99 

0.96 

0.91 

0.84 

0.76 

0.67 

0.58 

0.49 

0.42 

0.40 

0.5 

1.20 

1.19 

1.14 

1.07 

0.98 

0.86 

0.72 

0.57 

0.42 

0.27 

0.20 

a 

-  0.4 

r 

-  0.1 

0.70 

0.70 

0.68 

0.66 

0.64 

0.61 

0.58 

0.55 

0.52 

0.51 

0.50 

0.3 

0.90 

0.89 

0.86 

0.81 

0.75 

0.67 

0.58 

0.49 

0.40 

0.33 

0.30 

0.5 

1.10 

1.09 

1.05 

0.98 

0.89 

0.78 

0.65 

0.51 

0.35 

0.20 

0.10 

3A 
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TABLE  IV  (Continued) 


a 

«  0.5 

0° 

18° 

36° 

54° 

72° 

VO 

o 

o 

108° 

126° 

144° 

162° 

180° 

r 

»  0.1 

0.60 

0.60 

0.58 

0.56 

0.54 

0.51 

0.48 

0.45 

0.42 

0.41 

0.40 

0.3 

0.80 

0.79 

0.76 

0.72 

0.66 

0.58 

0.50 

0.40 

0.31 

0.23 

0.20 

0.5 

1.00 

0.99 

0.95 

0.89 

0.81 

0.71 

0.59 

0.45 

0.31 

0.16 

0.00 

a 

■  0.6 

r 

»  0.1 

0.50 

0.50 

0.48 

0.47 

0.44 

0.41 

0.38 

0.35 

0.32 

0.31 

0.30 

0.3 

0.70 

0.69 

0.67 

0.63 

0.60 

0.50 

0.42 

0.33 

0.24 

0.15 

0.10 

0.5 

0.90 

0.89 

0.86 

0.80 

0.73 

0.64 

0.54 

0.42 

0.29 

0.17 

0.10 

a 

=  0.7 

r 

»  0.1 

0.40 

0.40 

0.39 

0.38 

0.34 

0.32 

0.29 

0.26 

0.23 

0.21 

0.20 

0.3 

0.  oO 

0.59 

0.57 

0.53 

0.49 

C.42 

0.35 

0.27 

0.19 

0.09 

0.00 

0.5 

0.80 

0.79 

0.76 

0.72 

0.66 

0.58 

0.50 

0.40 

0.31 

0.23 

0.20 

a 

=  0.8 

r 

=  0.1 

0.30 

0.30 

0.29 

0.27 

0.25 

0.22 

0.19 

0.16 

0.13 

0.11 

0.10 

0.3 

0.50 

0.49 

0.48 

0.45 

0.41 

0.36 

0.30 

0.24 

0.18 

0.13 

0.10 

0.5 

0.70 

0.69 

0.67 

0.64 

0.59 

0.54 

0.48 

0.42 

0.36 

0.32 

0.30 

a 

-  0.9 

r 

-  0.1 

0.20 

0.20 

0.19 

0.18 

0.16 

0.14 

0.12 

0.09 

0.06 

0.03 

0.00 

0.3 

0.40 

0.40 

0.39 

0.38 

0.34 

0.32 

0.29 

0.25 

0.23 

0.21 

0.20 

0.5 

0.60 

0.60 

0.58 

0.56 

0.54 

0.51 

0.48 

0.45 

0.42 

0.41 

0.40 
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